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ABSTRACT 

A new open-loop heat flux control technique has been developed to conduct transient thermal testing of thick, thermally- 
conductive aerospace structures. This technique uses calibration of the radiant heater system power level as a function of heat 
flux, predicted aerodynamic heat flux, and the properties of an instrumented test article. An iterative process was used to 
generate open-loop heater power profiles prior to each transient thermal test. Differences between the measured and predicted 
surface temperatures were used to refine the heater power level command profiles through the iteration process. This iteration 
process has reduced the effects of environmental and test system design factors, which are normally compensated for by 
closed-loop temperature control, to acceptable levels. The final revised heater power profiles resulted in measured temperature 
time histories which deviated less than 25 °F from the predicted surface temperatures. 
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1. INTRODUCTION 


Ground thermal testing is often required to establish the flight worthiness of, or conduct research on structural components 

flight conditions generate thermal gradients in the structure, which only transient tests can accurately simulate. 

Transient thermal tests have been conducted on a variety of aerospace vehicles, or their components, including the X- 15 
horizontal stabilizer, 1 YF-12 aircraft, 2 and Space Shuttle elevon seals. 3 These structures had thin skins or were therma y 
lated and a traditional closed-loop temperature control system was used to control surface temperature. Temperatures on 
test articles were generally controlled to within ±15 °F of the desired temperature profile. 

open-loop heat flux control technique was developed for the transient thermal test of the PHYSX wing glove. 

This paper describes the open-loop heat flux control technique that was developed for use on the leading ; edge 
wing gk>r Background information is provided on the facility, test systems, and test article used du , nn g*e PHYSX tes ^ 
brief discussion of traditional feedback temperature control, as applied to transient thermal testing, an rep 
also^Dresemed 11 Experimental results obtained using the open-loop heat flux control method are discussed. The paper close 

improvements of the open-loop heat flux control method. Note that use of trade nam o 
names of manufacturers in this document does not constitute an official endorsement of such products or manufacturers, eithe 
expressed or implied, by the National Aeronautics and Space Administration. 

2. BACKGROUND 

The Flight Loads Laboratory (FLL, formerly known as the Flight Loads Research Facility, 5 and Thermostmctures Research 
Facility ^located at the National Aeronautics and Space Administration (NASA), Dryden Right Research Center, Edwards, 

designed to perform combined thermal 2 and mechanical loading of aerospace structures. 


®Pegasus is a registered trademark of Orbital Sciences Corporation, Fairfax, Virginia. 
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This section provides background information on: (1) the systems (including a thermal control system and radiant heater 
systems) used in the FLL for transient thermal testing of aerospace structures, (2) how transient temperature profiles are 
obtained, and (3) the PHYSX test article which sparked the development of the open-loop heat flux control method. 

2.1 Thermal control 

Temperature feedback control is used in many research and industrial applications. A wide variety of controllers are available 
with set point or multiple segment ramp and soak temperature control. However, the data acquisition and control system 
(DACS) in the FLL was specifically designed to perform a dynamic thermal simulation of a high-speed maneuvering aerospace 
vehicle. Figure 1 shows a schematic of this system, which was tailored to test thin wall structures, such as aircraft wing and 
fuselage sections. Temperature feedback control is implemented through an adaptive proportional integral algorithm. 8 Critical 
test article and heater dependent control parameters for the algorithm are obtained from a specimen identification test 
performed prior to the transient thermal test. 

The FLL DACS is capable of controlling up to 512 zones. Each zone consists of a portion of the radiant heater which illumi- 
nates a certain area of the test article surface. A thermocouple is installed on the test article within each zone and is used as 
feedback to the control system. The DACS forces the temperature in each zone to follow a prescribed temperature time history 
(or profile). Each zone is assigned a profile which may change slope as often as once per second. Open-loop power commands 
may also be prescribed for all or part of the profile of a zone. 

The electric power supplied to the radiant heaters is regulated by zero-crossover power controllers. This type of controller 
turns each cycle of 60-Hz alternating current on or off at the moment between cycles when the line voltage crosses zero (fig. 2). 
This control method minimizes electrical noise but causes pulsed heating that is especially evident at low power levels. The 
voltage and current supplied to the heater during each on cycle is set by the transformer which is feeding the power controller 
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Figure 1 Feedback temperature control system in the Flight Loads Laboratory. 
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correspond directly to the number of power cycles turned on in each 0.25-sec update interva . 

2.2 Radiant heaters 

to as quartz lamps, are commercially available from major lamp manufacturers .n a variety of wattages and lengths. 

The reflectors are used to contain the radiant energy within the test system and direct as much energy as possible to the test 
article surface These reflectors are manufactured from either fused silica ceramics or metals. Uncooled ceramic reflectors ca 
operate for long periods of time at temperatures in excess of 2000 °F. Metallic reflectors (usually aluminum °^ s ^n less f ee ) 
are used uncooled for testing at low temperatures «500 °F) and low heat flux «10 Btu/ft /sec). 
cooled reflectors, or air- and water-cooled reflectors are used in tests requinng test artic e empera u p 

types of reflectors are also used in tests requiring heat fluxes as high as 100 Btu/ft /sec. 


Radiation fences are used at zone boundaries to minimize 
typically made from fused silica ceramics. The ceramics 


the influence of a zone on adjacent zones (i.e., crosstalk) and are 
offer several benefits including tolerance of temperatures up to 
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Figure 3. PHYSX wing glove radiant heating system, removed from test article. 


3000 °F and minimal thermal expansion. Radiation fences are particularly critical when the heat fluxes being applied by 
adjacent zones are significantly different. The heater for the PHYSX wing glove (fig. 3) was designed with very narrow zones, 
approximately 1 in. wide, at the leading edge (left side of the figure) because the heat flux changes significantly from 
27 Btu/ft 2 /sec at the leading edge (left-most lamps in the figure) to 6 Btu/ft 2 /sec two zones (about 4 in.) away. 

2.3 Transient temperature profile generation 

Transient temperature profiles may be obtained from either direct flight measurements or through an aerodynamic heating 
analysis. To begin the analysis, known flight time histories of altitude, Mach number, and angle of attack for the vehicle are 
used as input parameters for an aerodynamic heating program. The program then calculates transient time histories of surface 
temperatures, heating rates, heat transfer coefficients, skin friction, and surface static pressures for selected locations on the 
vehicle structure. 
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The transient aerodynamic heating program used at NASA Dryden is an in-house program called THEOSKIN. Each location 
on the vehicle structure is modeled by a single thermal mass, or lumped capacitance element. Hus model assumes a urnform 
temperature throughout the element although the actual structure may have a nonuniform temperature distnbution. The thi 
ness P of the lumped capacitance element must be adjusted so that the element temperature in THEOSKIN represents the surface 

temperature of the real structure. 

Finite-difference thermal models at various locations of interest on the vehicle structure are then developed using the 
Lockheed Thermal Analyzer 9 (LTA). The heating rates calculated by THEOSKIN are applied to the LTAthermid models that 
are used to calculate surface and structural temperature time histories. Peak surface temperature values obtained from LTA are 
compared to the peak temperature obtained from THEOSKIN for each location. If the temperature deviations are greater than 
a few degrees, the element thicknesses used in the THEOSKIN calculations are adjusted. After the element thicknesses are 
adjusted and new heating rates are calculated, the thermal models are executed again with the newer heating ra^^d Ae peak 
temperatures compared. This iterative process continues until the peak surface temperatures calculated in THEOSKIN are 
within a few degrees of the peak LTA surface temperature at each location. The predicted surface temperature time histones 
from LTA are then used as the required surface temperature control profiles. 


2*4 PHYSX wing glove test article 

The PHYSX wing glove 4 is a flight test article designed to gather hypersonic boundary-layer transition data at speeds up 
to Mach 8. The flight test article is mounted on the right wing of a Pegasus air launched space booster (fig. 4). The glove is 



Figure 4. Plan view of Pegasus® booster and PHYSX wing glove. 
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manufactured from structural steel with the upper and lower surfaces being thin skins (=0.1 in. thick) and the leading edge 
being a thick (nominally 0.5 in.) structure in order to act as a heat sink. Numerous attachment fittings were machined into the 
leading edge, which more than doubled the local thickness of the part. The leading edge is rigidly attached to the upper and 
lower surface thin skins. The glove must maintain a prescribed shape during the Pegasus® first stage flight, beginning at an 
altitude of 38,000 ft and Mach 0.8 and progressing to an altitude of 200,000 ft and Mach 8 only 80 sec later. The launch 
temperature of the glove is -30 °F. Temperatures on the leading edge reach 500 °F by the end of the 80-sec flight. The peak heat 
flux predicted on the glove ranges from 27 Btu/ft 2 /sec at the leading edge to 6 Btu/ft 2 /sec 4 in. behind the leading edge. The 
peak predicted heat flux at the trailing edge of the upper surface (48 in. behind the leading edge) is 3 Btu/ft 2 /sec. 

A transient thermal ground test was required to validate the PHYSX wing glove design when exposed to the predicted aero- 
dynamic heating environment. The thermal gradients on the surface and through the thickness of the leading edge are the 
primary causes of thermal stress in the wing glove. Therefore, accurate, real-time simulation of the leading-edge surface 
temperature time histories were required during the transient thermal test. A duplicate wing glove dedicated to the transient 
thermal test was manufactured and attached to a mock-up of the Pegasus® wing. The glove surface temperatures predicted by 
an aerodynamic heating analysis of the Pegasus® vehicle were to be imposed on the glove by the F LL radiant heating systems 
using closed-loop temperature control. Instrumentation on the test article included deflection potentiometers, strain gages, 
thermocouples, and heat flux gages. Deflection potentiometers, strain gages, and thermocouples were placed at various 
locations on the structure to obtain data for comparison to structural analyses. Thirty-nine thermocouples were placed on the 
exterior surface of the wing glove for use in feedback temperature control. 

The heater for the PHYSX glove thermal test (fig. 3) was divided into 39 zones. The zones were grouped into rows with each 
row containing one to four zones, depending on the test article geometry and heating requirements. Each leading edge row had 
three roughly equal area zones. Figure 3 shows that the zones near the leading edge were smaller and had denser lamp spacing 
in order to apply the wide range of heating rates predicted for the leading edge. 

A typical test began by precooling the test article to -30 °F to simulate the cold soak condition at Pegasus® launch. This 
process generally took 45 minutes to 1 hour to complete, and resulted in significant frost buildup on the test article. Once the 
precooling process was complete and the glove attained a uniform -30 °F temperature, the 80-sec transient thermal test 
was executed. 

The FLL closed-loop temperature control system performed well on the thin skin areas of the glove (fig. 5).The control 
system accurately compensated for the presence of the frost buildup by increasing the power to the heaters as the frost removed 
heat from the test article by melting and boiling away. The closed-loop temperature control system generated temperatures 
within ±10 °F of the desired profiles on the thin skins. 

The closed-loop temperature control system did not provide adequate control in three rows which heated the leading edge 
structure because of the leading edge thickness and the heating rates being applied. These rows included 9 of the 39 control 
zones used in the test. The controlled temperature in the nine leading edge zones began to oscillate as the heating rate 
approached 5 °F/sec (fig. 6). As the heating rate exceeded 5 °F/sec, the oscillations increased to unacceptable levels. The devi- 
ations from the profile eventually exceeded limits, which caused an automatic shutdown of the heater system. This oscillating, 
diverging behavior was not an acceptable representation of the predicted flight profile. Because closed-loop temperature 
control was successful on the remaining 30 zones, this method was used to control them for the duration of the test program. 

3. OPEN-LOOP HEAT FLUX CONTROL 

An alternative to closed-loop temperature control for the leading edge upper, lower, and center rows of the heater is to direct- 
ly control predicted absorbed heat flux at the surface of the test article. Unfortunately, the DACS thermal control software was 
optimized for closed-loop temperature control of thin skin test articles and could not be easily modified to use a closed-loop 
heat flux control algorithm. 

Open-loop programming of the heater power levels was the only near-term alternative to temperature feedback control on the 
DACS. This would allow delivery of a prescribed incident heat flux profile but would not provide real-time compensation of 
changes in the absorbed heat flux caused by environmental factors. Therefore, the open-loop heat flux control method is neces- 
sarily an iterative process. It requires knowledge of the predicted and measured surface temperature profiles {T a (t) and T m (t ) ), 
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Figure 5. Temperature control of a thin skin zone. 



Figure 6. Temperature control of a leading edge zone. 


the relationship between the heater power commanded and the absorbed heat flux obtained (represented by ,P q ), the predicted 
^“flux profiles (, (,) ), and knowledge of the thennal properues of the test article, such as the densuy (p). spectfic 


heat (C p ), and surface thickness (t). 
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3.1 Pretest activities 


Figure 7 shows a flow chart summarizing the open-loop heat flux control technique described in this section. The process 
begins with the installation of calibrated heat flux gages on the surface of the test article. One heat flux gage was installed on 
the glove surface at the spanwise center of the leading edge upper, lower, and center rows (fig. 8) to directly measure the heat 
flux generated by the quartz-lamp heaters. It was assumed that the heat flux delivered at a given power level was constant along 
each row although each row consisted of three control zones. Therefore, only one heat flux gage was installed in the center 
zone of each row. Thermopile heat flux gages were installed in the leading edge upper and lower rows. A Gardon heat flux 
gage, capable of operating at temperatures in excess of 500 °F, was used in the leading edge center row. At this location, 
expected temperatures exceeded the capabilities of the thermopile heat flux gages available for use on the PHYSX wing 
glove test. 

To apply the proper heat flux to the test article surface, the relationship between the heater power level commanded by the 
DACS and the heat flux actually produced by the heater was required. This relationship was assumed to be linear, based on past 


Install calibrated heat flux gages 


I 


Perform heater calibration test 


p _ heater power 
*1 ~ heat flux 


Program heater power, P(t) 


P(t) = P q q(t) 


Start 


Perform aero/thermal analysis 


• Determine proper heat flux 

time history, q(t) 

• Predict glove temperature 

response, T a (t) 


Perform transient heating test 


• Measure glove temperature 
time history, T m (t) 


* 



Update heat flux 
until convergence 

q(t) = q(t) + (Aq) n (t) 



Figure 7. Open-loop heat flux control method. 
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Figure 8. Heat flux gages on PHYSX wing glove. 



970707 


experience in the FLL. The validity of this assumption is addressed in section 4.3. To obtain this relationship, a heater 
calibration test was performed at a constant power level. The level chosen was high enough to minimize the sawtooth heat flux 
profile caused by zero-crossover power control, yet low enough to limit temperatures and stresses in the glove to acceptable 
levels Power level 5 met this requirement and was entered into the control computer. Each row of the heater was fired individ- 
ually to obtain its power level per heat flux calibration (P q ). Table 1 shows the calibration values obtained during the test. 
Crosstalk between the rows was less than 1 Btu/ft^/sec at power level 5 (fig. 9). 


Table 1 . Heater calibrations. 


- 

Heater calibration 

Leading edge 

power level 
2 

heater row 

Btu/ft /sec 

lower 

0.87 

center 

0.4 

upper 

0.91 


The calibration test revealed significant temperature dependence in the output of the Gardon heat flux gage in ^the leading 
edge center row (fig. 9). The temperature dependence was observed as the heat flux gage temperature increased during e 
10-sec application of a constant power level. The indicated heat flux dropped from 14 Btu/ft /sec to 10.3 Btu/ft /sec, although 
the power level (i.e. heat flux) was constant throughout the calibration test. Thermopile heat flux gages used in the upper and 
lower rows indicated a constant 1 Btu/ft 2 /sec while the center row was firing. The thermopile heat flux gages produced essen- 
tially constant heat flux measurements while their own rows fired as well. This confirms that the center row lamps emitte a 
constant heat flux and the indicated reduction in heat flux was a characteristic of the Gardon heat flux gage. 


Another key step in the open-loop heat flux method is determining the heating rates required to accurately simulate the 
thermal environment. Aerodynamic heating rates calculated for each zone from the THEOSKIN program were used to deter- 
mine the net heat flux profiles (q(t ) ) required for open-loop heat flux control. Because these are the same heating rates used in 
determining the transient surface temperature time histories discussed in section 2.3, there was little additional effort required 
to generate these heat flux profiles. 
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Time, sec 

Figure 9. Measured heat flux during leading edge center row calibration test, center zone. 
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The heater power commanded by the DA CS during the calibration test was divided by the actual heat flux obtained at the 
surface to obtain the heater calibration coefficient (P ) necessary for the control method. The heater power time histories were 
then generated by multiplying the predicted heat flux profile for a particular row by the heater calibration of that row. 

F(r) = P q q{t) y d) 

where P(t) is the power level programmed into the DACS, as a function of time. The crosstalk between the rows was not taken 
into account while generating power level profiles, because it was a relatively small component of the total heat flux being 
applied. An iteration process was used to account for all crosstalk as well as other environmental effects that were known to be 
present but were unpredictable. 

3.2 Iteration process 

After determining the necessary commands to the heaters, an initial test was performed and the transient temperature distri- 
butions at each zone on the glove were measured. Following the test, these measured temperatures (T (/) ), were compared to 
what was analytically predicted <J a (t ) ) for each second, and a differential heat flux profile for that time interval was calculated 
using the following equation (appendix A) 


A<7„(0 = PC i 




At 


At 


(2) 


where A q n (t) = differential heat flux at iterationn^expressed ^as a function of profile time; p and C p represent the density and 
the specific heat of the test article, respectively; ~ and — — are the predicted and measured temperature rise rates at 
the heated surface of the test article, respectively; and x is the test article thickness. Note that equation (2) is based on a uniform 
temperature change per unit thickness. With thick structures subjected to transient heating, the temperature distribution through 
the thickness may not be constant. The full thickness of the test article at the leading edge center row (x = 0.5 in.) was originally 
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»«d in equation (2). An effective thickness (,*) was then entphicall, detetnhned based on a second test in order to increase the 
me of convergence. The differential heat flux profile in equation (2) was then used to update heat flux profile, ,„<r). 




(3) 


This process was repeated until die ^ /^^Rt^lySTSc^t SIC 
temperature profiles within the bounds of «%£££% J^ff“SS«ts pmanteter (x*> can be analyfically 
not used in this test program, appendix B pr ,-mnerature can be acquired during the test. Determining an 

^t;b^o!x* n ur»gTe' 'IlyMnteld proposed in appendix B may increase die convergence rate in the item 

lion process. 

3.3 Limitations 

The open-loop heat flux “ "ccomTfoTtecTu^of ^ tab dtdS'^Sons. Wind and temperature affected 

t^PHY^ g^ve tm^becairw^oors a "d pppyg^'glov^lest adde^H^t^the^msient'hmlng test 

glove. The precooling process also caused m bu.W up on the PHT« 8^^ ^P ^ ^ temperature as 

The frost accumulation increased the thermal t effects on the PHYSX glove were relatively small 

the melted frost boiled off the glove. The wmd and ambient temperature C too could not be determined 
compared to the effect of frost buildup. However, the ^nT^Sve^s us^ the same tempera- 

because dtese facom could no. be independently venal. ™ IT.™ changes The measured sur- 

ture and power level profiles in order to quanu y e ove P as 25 °F when the results of the two tests were 

— ~ • - — — - 

weeks, months, or longer, unless environmental factors can be well controlled. 

Tire FLL test system itself presents a set of ^-a^sov^wer c^Tres^ flux profile 

system (both facility and test specific) design. In the FLL, zero . crosso p which 

at low power levels and discretized power command. The J calibra - 

PHYSX wing glove, difficult, 
effect to acceptable levels. 

The test article may also impose limits on the method. If the transient thermal test has a significant chance of damaging the 
test article or is planned as a test to i failure, the^ability^to of dam- 

aXTTesr^icle duringTaf time. A real-time modification of the open-loop heat flux control power level profiles or 
improved temperature feedback control would greatly reduce delations from the desired profile during such one-shot tests. 

Th, final difficulty with the method is heat flux sensor calibration uncertainty. One should expect large uncertainties in heat 
flux^^ure^nte^iwnflie cu^i^state of the art. The iteration process will eventually correct for this uncertainty, but a more 
accurate heat flux measurement leads to fewer iterations to converge on the desired temperature profile. 

4, EXPERIMENTAL RESULTS AND DISCUSSION 

Test data obtained on the open-loop heat flux control method can be grouped into three parts^ First, data from a zon ewhere 
a heat flux gage was installed will be discussed. Second, data from a zone that did not have a heat flux gage are presente 
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highlight the ability of the iterative method to compensate for heater calibration errors. Finally, data are presented to validate 
the single power level calibration process. 

4.1 Zone with heat flux gage installed 

Test data from a zone with a heat flux gage installed is shown in figure 10. Note that tests performed to support the iteration 
process were intentionally stopped at 60 sec to prevent possible damage to the test article. Once the iteration process converged 
in the first 60 sec, the full 80-sec test was run with no iteration performed between 60 and 80 sec. All tests began with the zone 
under temperature feedback control and switched to open-loop control at 12 sec. This allowed closed-loop temperature control 
to compensate for environmental effects at the beginning of the profile and provided a repeatable starting point for the open- 
loop portion of each profile. Temperatures within ±25 °F of the desired profile were generated by the test performed after iter- 
ation, which is within the expected day-to-day variation caused by environmental factors. The fact that the measured tempera- 
ture remained reasonably close to the predicted temperature profile after the 60-sec point indicates that the heater calibration 
used for this zone was accurate. In addition, environmental factors, such as frost, were the cause of the first test deviation from 
the desired profile. 

4.2 Zone without heat flux gage installed 

The iteration process was particularly valuable in the open-loop controlled zones without heat flux gages. Errors in the 
assumed heater calibration can be large, because heat flux was not measured in these zones during the heater calibration test. 
Figure 1 1 shows the results obtained during transient testing. The measured temperatures were within the error band that was 
expected because of environmental effects. The deviation from the desired temperature grew after the 60-sec point when no 
iteration was possible. The open-loop heat flux control method provided an acceptable temperature profile even in zones where 
heat flux could not be directly measured. 

4.3 Heater calibration validation 

The heater calibration calculation was performed for the leading edge center row using heat flux and power level data 
measured during the final transient thermal test (fig. 12). The original calibration was performed at a constant power level 5. 
Significant scatter and nonlinearity were noted in the data collected during the transient test below power level 5. Much of the 
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Figure 10. Results of open-loop process, zone with heat flux gage. 
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Figure 12. Heater calibration calculated from transient test data. 


nonlinearity was caused by variation of the indicated heat flux with heat flux gage temperature (section 4.1). The scatter was 
caused primarily by the zero-crossover power control. The calibration calculation assumes that a constant heat flux value is 
generated at each power level, not the sawtooth pattern created by low power level zero-crossover control. A wide range of 
calibration values at low power levels results from the sawtooth heat flux pattern. 
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The transient heater calibration data taken above power level 5 shows that the linear assumption for the heater calibration is 
reasonable. The calibration in this power level range is nominally 0.5 power level/Btu/ft 2 /sec, somewhat above the calibration 
shown in table 1. This offset is due to the temperature-dependent output of the Gardon heat flux gage. 

5. FUTURE WORK AND RECOMMENDATIONS 

The method developed in this paper produced transient temperatures on the glove that were well within the program require- 
ments. However, much of what was developed for the PHYSX represented a work-around solution given the time and budget 
constraints of the project, as well as limitations with the current data acquisition and control system used in the tests. The open- 
loop heat flux control technique was developed primarily because of a limitation with the current thermal control system in 
controlling the temperature of thick structures. Therefore, the first step for future work is to implement some of the lessons 
learned from this work into a new DACS. A new system with higher feedback rates in the control loop may improve the viabil- 
ity of temperature control of thick structures. Design and development of such a system is underway. 

In many applications, temperature feedback control is not a viable option and the actual heat flux environment must be sim- 
ulated. To this end, steps are currently being taken to program the algorithm developed in this paper into the next generation 
DACS so that adjustments to the heat flux profile can be made real time, instead of between tests. To avoid some of the limita- 
tions noted earlier with the present algorithm, other options are being explored to obtain accurate calculations of the differential 
heat flux. 


Front and back surface temperature measurements should be made whenever possible during transient tests of thick 
structures. With such measurements, it is possible to integrate a one-dimensional finite-difference model in the iteration pro- 
cess. With this approach, the front and back surface temperatures measured during the test would be prescribed at boundary 
nodes of the model, and the nonlinear distribution through thickness of the surface would be calculated as well as an estimate 
of the absorbed surface heat flux. The differential heat flux required would be determined by subtracting the estimated surface 
heat flux from the absorbed heat flux predicted by the aerothermal analysis. With recent advancement in real-time processing, 
this process could be implemented real time instead of having to rely on posttest processing. 

It will not always be possible to measure the back-side temperature of the test structure. In such a case, one must either 
assume an appropriate boundary condition at the back side of the structure in order to determine the absorbed heat flux analyt- 
ically or make the constant temperature assumption. Reasonable values for T* for use in equation (2) can be determined from a 
pretest thermal analysis as shown in appendix B. 

6. CONCLUDING REMARKS 

The new open-loop heat flux control technique that has been developed to conduct transient thermal testing of thick, 
thermally conductive aerospace structures uses calibration of the radiant heater system power level as a function of heat flux, 
predicted aerodynamic heat flux, and the properties of an instrumented test article. Open-loop heat flux control provides a 
reasonable alternative when the test article, test conditions, or both cause unstable temperature feedback control during tran- 
sient thermal tests. Data from the Pegasus® Hypersonic Experiment (PHYSX) wing glove thermal test show that the method is 
successful in a challenging thermal test environment. Heat flux time histories and the heat flux per power level heater calibra- 
tion must be determined prior to transient testing, in addition to the information usually required for temperature feedback 
control. Steady-state calibrations of the radiant heating system can be performed with sufficient accuracy to begin the iteration 
process. The iterative modification of the power level (i.e. heat flux) profiles is critical in order to account for various effects 
introduced by the environment, test systems, test article, and instrumentation. 

Heat flux gages need not be installed in every zone if multiple zones can be expected to deliver similar heat flux to the test 
article. The assumptions made during the heater calibration test were confirmed by transient test data. The open-loop heat flux 
control method generated transient surface temperature histories within ±25 °F of the desired profile on a thick, complex part 
exposed to large changes in applied heat flux over short distances. The iteration process was successful although a heat flux 
gage was installed in only the center zone of each spanwise row. The success of the open-loop heat flux control on the PHYSX 

wing glove warrants development of a real-time algorithm to provide increased accuracy during one-shot transient 
thermal tests. 
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7. APPENDIX A 


, , . . , flliv ( nprpccarv to ensure that the measured surface temperatures 

predicted by the analysis. This appendix develops the mathematical basis for 

this equation. 

When a structure of density (p), volume (V), and specific heat (C„) is heated transiently such that the temperature increases 
by'adifferential change in temperature (dT), the change in internal energy of the structure (du) is given as 

du = pVC p dT(x,t) (A1) 

to thickness T,'toe change in interna, energy of a one-dimenstonal eienren, wito uni, area in: 


t T 

u" = jdu"'dx = jp C p dTdx 
0 0 


(A2) 


where in" represents toe change in internal energy per uni, volume. If aU toe terms on toe right side of equation A2 are 
constant with respect to x, equation (A2) becomes 


u" = pC p xdT 


(A3) 


The heat transfer rate per uni, area, or hea, flux (,(0) can he determined simply by differentiating equation (A3) with respect 
to time. 


q{t) - P c p ‘ z dt 


(A4) 


Eouation (A4) can now be employed to determine the surface heat flux as a function of time. This heat flux profile must be 
applied to produce a temperature time history (7X0) on the surface of the structure, provided that the temperature distn u 
through the thickness of the structure remains constant. 

Figure A1 shows toe case encountered with toe initial testing of toe PHYSX glove in which the towntod I transient temper, 
ure at the heated surface (Tlx = 0, /)) differed from a desired surface profile calculated analytically (TJx - 0, 0). 


atUTe 

nomenclature in figure A1 is simplified with the following expressions: 

T(t) = T(x = 0,f) 

A T {S+l) = T(x = 0,t= 5 +l)-r(x = 0.r = s) 
Ar (5 + 1) = r(s+l)-f(0 


(A5a) 

(A5b) 

(A5c) 


From figure Al, the net change in heat flux (A q) required in the time interval between s and 5 + 1, (Ar (i + 1 } ) to equate the 


two curves is 


Ad At 


+ m _ 


p V 


/ A7' (i+1) AT (I+1)> ' 

a m 


At 


(i + 1 ) 


At 


.(r+D 


(A6) 
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8. APPENDIX B 

For many applications such as the heating of thin skin structures, the temperature drop through the skin is negligible, and 
equation (2) is appropriate. However, for an application in which a thick structure is being transiently heated, the temperature 
through the structure may not be constant; and therefore, this equation must be modified. The nature of the temperature varia- 
tion through the structure depends on the material properties, structural thickness, and transient heat flux profile applied. This 
appendix describes how the thickness parameter in equation (2) can be modified to help improve the convergence rate in the 
iteration process. This method also provides the basis for a closed-loop heat flux control method which is the subject of 
future work. 

In appendix A, the change in internal energy corresponding to a temperature change (dT) in a structure of specific heat ( C ), 
density (p), and thickness (x) is given in equation (A3). 

If the temperature change in equation (A3) is taken between the desired analytically determined surface temperature, 
(T a (x = 0)) and the measured surface temperature (T m (x = 0)) both taken between time 5 , and time s + 1, equation (A3) 
becomes: 


= p V 


r(*+l) 


£ +1) ) 


(Bl) 


where T 


<J+1) 


, and T 


(s+ 1) 


are both assumed constant in x (fig. B 1). 



Figure B 1 . 
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For the transient heating condition where the temperature distributions are not constant in x, as shown in figure B2, the 
change in internal energy in equation (A2) becomes: 


rx 


«' = pc, 




,(s+ 1) 


w-4 +l) w]^ 


LO 


(B2) 


Because equation (Bl) is a specific case of equation (B2), the change in internal energy determined from the two equations 
will generally be different (fig. B2). The goal is to define an effective thickness parameter (x*) which, when multiplied by the 
surface temperature, will produce the same area as the shaded region in figure B2. 



0 T X 

Thickness, In. ^ 

Figure B2. 


Redefining T in equation (B 1) as the effective thickness, (x*), setting equations (Bl) and (B2) equal, and solving for x*, 

T 


J(r<” 1) w-T< m * +1) w )* 

0 

(J + l)_ T (s + 1) 

1 a m 


(B3) 


which is shown graphically in figure B3. 



0 


T* 

Thickness, in. 

Figure B3. 
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Proper selection of T* will improve convergence. Two methods are suggested to select an effective thickness. The first case 
is when both the front and back surface temperatures (7^(0, t) and TJX, t)) of the heated structure are measured during the 
transient thermal test. The second case is when only the front surface temperature of the structure is known. For the first case, 
T^ S+1 \x) can be determined at discrete nodal locations through the structure thickness using an aerothermal analysis, as 
described in section 2.3. If the front and back surface temperatures (7*^(0, /) and T m ( t, f)) of the heated structure are mea- 
sured during the transient thermal test, then these measurements can be prescribed as boundary conditions in an analysis, the 
interior temperature distribution, (7^ + !) (x)) can be determined analytically, and equation (B3) can be solved. Therefore, an 
appropriate thickness parameter is determined which ensures that the change in internal energy of the isothermally heated 
effective structure is equivalent to the change in internal energy produced by transient temperature distribution in the 
actual structure. 

For the case where only the front surface temperature of the structures is known, further assumptions are in order. Assuming 
that the error between the desired and measured temperature profiles is caused by a greater uncertainty in the heat flux applied 
than in the material properties, and if the temperature distribution calculated from the test is reasonably close to the ^predicted 
values, then X* can be approximated strictly based on the analytically determined temperatures. Replacing T m (x) with 
T^\x) in equation (B3) yields: 


T* 




0 

(r+1) (*) 

1 a 1 a 


(B4) 
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